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Abstract. Here we prove a new lower bound on the indirect Coulomb energy in quantum 
mechanics, in terms of the single particle density of the system. The new universal lower 
bound is an alternative to the classical Licb Oxford bound (with a smaller constant, C = 
1.45 < Clo — 1-68) but involving an additive kinetic energy term of the single particle density 
as well. 



1. Introduction 

Consider a system of N particles with charges e\, . . . , > 0. In non relativistic quantum 
mechanics, this system is described by a normalized wavefunction, 

il>(xx,...,x N ; <?!,..., cr N ), (1) 

where x\, . . . , xn denote the coordinates of the particles (xi 6l 3 ), and <7j represent possible 
discrete quantum numbers, such as spin. The corresponding probability density function is 
given by, 

f(x 1 ,...,x N ) = ^ \ip(x u ...,x N ;a u ...,a N )\ 2 . (2) 

Cl,...,<TjV 

We define the charge density of particle i as, 

pi(x) = e % J f(x U X, X i+1 , . . . , X N )dXi, (3) 

where dxi means integration in all the particle coordinates, except the i-th. Finally, we then 
define the single particle density as, 

N 

p(x) = ^ Pi (x), (4) 
i=i 

which is a charge density for the whole system. Here, we assume a Coulomb interaction 
between the particles given as usual by 



l<i<j<N 1 



The expectation value of this interaction when the system is described by the wavefunction i/j 
can be simply expressed as, 

£ f JSZLdX, (6) 

I 1. J | 



i<j 
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where X — (xi, . . . , x^) and dX = dx\ . . . dx^. This expression can be decomposed into two 
parts, 

2 J \x-y\ 

where the first term, which it is usually called the direct part of the Coulomb energy, represents 
a classical expression for the electrostatic energy of a charge distribution p. On the other hand, 
the remainder E is all that has been missed by treating the system as such a distribution p, 
and it is known as the indirect part. In this article we will be interested in finding a lower 
bound for E. 



In 1930, Dirac j| gave the first approximation for the indirect Coulomb energy in terms of 
the single particle density. In the standard case, when e; = e, the absolute value of the charge 
of the electron (for all i), using an argument with plane waves, he approximated E by 



E 



-c D e 



2/3 



^ 3 dx, 



where c D = (3/4)(3/tt) 1 / 3 « 0.7386 (see, e.g., @, p. 299). 

The first rigorous lower bound for E was obtained by E.H. Lieb in 1979 H, using the 
Hardy-Littlewood Maximal Function [|13|]. There he found that, 



E > -8.52 



N 



1 4/3 >| 3 / 4 

dx 



i=\ 



p(x) 4/3 d 



X 



1/4 



(9) 



The constant 8.52 was substantially improved by E.H. Lieb and S. Oxford in 1981 [0, although 
the form of the lower bound was slightly different, 



E > -1.6 




N 



^eiPi{x) 



4/3 ^ V 2 

dx 



1/2 



(10) 



However, in the more standard case, when e\ 
to a similar form, 



p(x) 4/3 dx 

e^v = e, the two inequalities are reduced 

E>-C^l^^ (11) 

where C is a positive constant. The best value for C is unknown, but Lieb and Oxford 
proved that it is larger or equal than 1.234. The constant they obtained in general, 1.68, was 
found by first assigning to each particle a spherically symmetric charge distribution; however, 
they were not able to obtain the distribution that completely minimized the constant. Then 
Chan and Handy, in 1999 |l|, found a better value, 1.636, by optimizing numerically this 
distribution of charge. It is this last constant, as far as we know, that is the smallest value for 
C that has been found up-to-date. During the last thirty years, after the work of Lieb and 
Oxford 0, there has been a special interest in quantum chemistry in constructing corrections 
to the Lieb-Oxford term involving the gradient of the single particle density. This interest 
arises with the expectation that states with a relatively small kinetic energy have a smaller 
indirect part (see, e.g., 



12] , |14| and references therein). 
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Here we will provide a lower bound that has a smaller value for C than any of the previously 
mentioned results, but at the price of obtaining another term that involves a gradient of a 
power of the single particle density. Our main result is the following: 

Theorem 1.1. For any normalized wave function ip(xi, . . . ,xn) and any < a < 1 we have 
the estimate 

E(ip) > -1.4508(1 - a)-^ 3 e 2 ^ J R3 p(x)^ 3 dx 
-0.2382(1 - a)-V6 e 5/6 ((^ ^^p)) 1 ' 2 ( f R3 p( x y/ 3 dx) 1/2 - e(1.1781/a)(V& \p\Jp) (12) 
where 

(VpMVp) ■= [ \7p(k)\ 2 \^k\dk = -L f [ 'V^-V^)' 2 ^, (13) 
where f(k) denotes the Fourier-transform 

f(k) = l e- 2mk -*f(x)dx . 



Remarks 1.2. i) For many physical states, the contribution of the last two terms in ( 17^ ) is 
small compared with the contribution of the first term. We illustrate this fact in the Appendix; 

ii) For the second equality in fifty see, e.g., ; Section 7.12, equation (4), p. 184; 

Hi) It was already noticed by Lieb and Oxford (see the remark after equation (26), p. 261 
on JTJj, that somehow for uniform densities the Lieb-Oxford constant should be 1.45 instead 
of 1.68; 



iv) In the same vein, J. P. Perdew |ll| , by employing results for a uniform electron gas in 
its low density limit showed that in the Lieb-Oxford bound one ought to have C > 1.43 (see 

also, 

2. Proof of ([12]) 

In this section we give the proof of the new lower bound on E. Through out this section 
we set e = 1. The proof of this theorem hinges on the following lemma, originally due to L. 
Onsager PU[ . We quote it from || for the case where all the charges are equal. 

Lemma 2.1 (Onsager's lemma). Consider N unit point charges located at the distinct points 
Xi, . . . , xn- For each 1 < i < N let p Xi be a nonnegative, bounded function that is radially 
symmetric about he point Xi and whose integral J p Xi {x)dx = 1. Then for any non-negative 
integrable function p we have the inequality 

. N N 

E i ._ I > - D (pi p) + 2 E D (p> p^ - E D ^ p^ ■ ( 14 ) 

%<] •> i=l i=l 

For the simple proof we refer the reader to ||. As in [0] for each 1 < i < N we set 

PxX x ) = X 3 Pi,{xi)p{Xp{xi) 1/3 (x - Xi)) (15) 

where p is a non-negative, bounded, radial function that integrates to one. Applying Onsager's 
Lemma with p = p^ to the indirect term 

E{ip) = <v> E i .1 - D (p^ p*) ( 16 ) 



i<j 
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one finds after simple calculations 



E(ip)>- f dydxp(y)Xp(x) 4/3 R(Xp(x) 1/3 (x-y))- XD(p,p) [ p{x) 4/3 dx (17) 

</IR 3 xR 3 </R 3 



Here 



where 



We set 



and 



R(t) = - t - 4>(t) (18) 
<P(t) = [ min ( j, -L ) p{y)dy . (19) 
F = [ dydxp{y)\p{xY lz R{\p{x) llz (x-y)), 

</R 3 xR 3 

k — [ dyix m^f^i, 

Jl3xl3 \ x — y\ 

L= f p(x) 4 / 3 dx, 

M i = ^T \z\ 2 p(z)dz, 

6 J\z\<l 

1/2 



M 2 = ( / \z\ 4 R{z) 2 dz 

With this notation we have: 
Lemma 2.2. 



" A 2 A 3 / 2 



F l/2 + ^2 L l/2 

A 



(20) 



Proof. We write 

J M 3 XM 3 ^rfxp(y)Ap(x) 4 / 3 J R(Ap(a;) 1 / 3 (a; - y)) = J R3xR3 dydxXp(x) 7 / 3 R(Xp(xy/ 3 (x - y)) 
+ | M3xM 3 dydx[p(y)^ - p(x)^][p(y)^ + p(x)^}Xp(x)^R(Xp(x)^(x - y)) (21) 
The first term on the right side can be computed and yields 

/ dydx\p{xj l l*R{\p{xyi\x-y)) = ^ / p{x) A ' 3 dx I R{\z\)dz (22) 

Now 

/ R(\z\)dz = 2ir - (f)(z)dz = 2n - <&(z)p(z)dz 

JR :i J\z\<l Jr 3 

where $ is the potential of a uniform charge supported in the unit ball with total charge 4n/3. 
This function can be readily computed to be 



and we get 



f + f(l-|^| 2 ) if |*| < 1 

4vr 1 



3 \z\ 



if \z\ > 1, 



/ R(\z\)dz = f \z\ 2 p(z)dz. 
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Thus, we have that 

f dydx\p(x) 7/3 R(\p(x) 1/3 (x~y)) = 4 f p(x) 4/3 dx^-[ \z\ 2 p,{z)dz 

The second term in (|21~D reads 

dydxipiy) 1 ' 2 - p(x) 1 />( Z/ ) 1 / 2 Ap(x) 4 / 3 J R(Ap(x) 1 / 3 (x - y)) 

<R 3 

dydx[p(y) 1/2 - p{x) 1/2 }p{x) 1/2 \p{xf 3 R{\p{x) l ' 3 (x - y)) 



(23) 



+ 



(24) 



Further, 



/ R3xR 3 dydx[p{yf/ 2 - p(a;) 1 />( y )V2 Ap(x) 4/3 jR(Ap(a;) i/3 (x _ y)) 
Jr3 xR3 rfj/fc W^ 1/2 V (?/)^|x - y| 2 Ap(x) 4 / 3 J R(Ap(x) 1 /3 (x _ y)) 



(25) 



and by Schwarz's inequality this is bounded above by 



f ^,^ r [pfa) 1/2 -pW 1/2 ] 2 



1/2 



X 



/ R3xK 3 dydxp{y)\ 2 p{xf' 3 \x - y\ 4 R{\p{x) 1 ' 3 {x - y)f 



1/2 



(26) 



Here we note that the function R(t) = for t > 1 and hence we can assume that the domain 
of integration is such that Ap(x) 1//3 |x — y\ < 1. Moreover, 



i2(f) < - . 



Thus, the second factor is bounded by 

1 



A3/2 



^xp(|/)Ap(x) 4/3 J R(Ap(x) 1/3 (x - y)) 



1/2 



i.e., 



< 



A3/2 



/ R3xR3 dydx[p{yf' 2 - P (xy/ 2 ]p( y y/ 2 Xp(xr/ 3 R(Xp(xy/ 3 (x - y)) 
Jr3 x r3 dydx ^lzffi 1 ^ 2 } [j M3xM 3 dydxp{y)\p{xYl 3 R{Xp{xfl\x - y))] V2 (27) 
We estimate the second term in (|24| ) in a similar fashion. Schwarz's inequality leads to 
/ K 3 XR3 dydx[p{yf/ 2 - p{xyi 2 ]p{xyi 2 \p{xyi 3 R{\p{xy/\x - y)) 



< 



Zr3 x r3 dydx 



\x-y\* 



1/2 



X 



J M 3 XR3 dydx\x - y\*\ 2 p(x) n / 3 R(\p(x)^(x - y)f 



1/2 



Changing variables y z = \p{x) x l 3 {x — y) yields for the second factor 



A -5/2 



p(x) 4/3 dx / \z\ A R{zfdz 



q 1/2 



(28) 

(29) 
□ 
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Proof of Theorem 1.1. For any a > we have that 



1 -K l ' 2 F 1 ' 2 <JL- + aF, 



A 3 / 2 - AaX 3 

and hence, 

Mi M 2 1/9 1/9 K 

-(l-a)A 2 (l-a)A 5 /2 4(1 - a)«A 3 ' 

Thus, the absolute value of the indirect term is bounded above by 



\E\ < 



r M 1 + A J D(/i, y u) 



(1 — ajA 5 / 2 4(1 — a)a\ 6 



(31) 



— a)A 

Now we optimize the first term with respect to A and obtain the upper bound 

3M 1 1/3 £( / i,p) 2 / 3 1 M 2 £>(/i, /i) 5 / 6 1/2 1/2 2. ( D(ji,n) \ 

22/3(1 _ a )l/3 + (1-«)V6 (2M!) 5 /6 4a V 2M i / 

If we optimize the expression 

0(Ai, /i) 2/3 Mi /3 = ^) 2 / 3 | ^ |*| V3 

with respect to the function fi we find that the optimizing \x is the uniform distribution on 
the unit ball. Simple calculations yield for this distribution 



and moreover M 2 = a/237t/2310. Thus, 

E>-± V3 L - „ 1 ,*F>J*j*. K V>LV> - . (32) 



10^3(1 -a)/ (l-a)V6 v 47r y V 2310 IQna 

□ 

3. Appendix 

In this Appendix we illustrate how for some states, the first term in (|T2| ) can be much 
larger than the second. In our first example we will consider a single particle density given 
by an exponential, normalized to iV (the number of particles). We show that in this case, 
the quotient between the first and the last term goes like 1/iV 1 / 3 , so that in the large particle 
number limit, the last two terms in flTjD are negligible. 

A function / : R 3 -> C is said to be in H l / 2 (B?) if / e L 2 (R 3 ) and \k\ 1/2 f(k) is in L 2 (M 3 ). 
For functions in this set, the following identity holds, 

= [f lfia i ) ~ f( ;! )l2 dxdy= f \2irk\\J{k)\ 2 dk, (33) 
\x — y\ J 



2tt 2 



(see, e.g., ||, Section 7.12, equation (4), p. 184). 
It suffices to prove that the quotient 

k(Vp) _ 

/ p 4 / 3 dx 



Q{p) (34) 
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can be made arbitrarily small. To this end, we consider the function, 

p(x) = ae-W. (35) 



Since f p(x) dx = N, 



Certainly yfp G L 2 (R 3 ), and 



so y/p G # 1/2 (M 3 ). Therefore, 



On the other hand, 



So we conclude that 



87roy /3 



& =(^) • (36) 



Tm = ^f h , (37) 

V V (|4ttA;| 2 + 6 2 ) 2 V ; 



4/ o , 27a 4 / 3 7r . . 

P 4/8 cfa = - § p-. (39) 



W = 3 %2 W i/3 » ( 4 °) 



and taking the limit iV — > oo, the desired result is obtained. 

As a second example, which is more relevant from the physical point of view, we use the 
Thomas-Fermi density of an atom, p^p. It is well known that the Thomas-Fermi density of 
a neutral atom of nuclear charge Z satisfies the following scaling, 

p TF (x) = Z 2 f(Z^x), (41) 

where / denotes the Thomas-Fermi density in the case Z — 1. Using this simple scaling 
relation, one can immediately check that 



J p 4 TF /3dx = cZ 5/3 , 



for some positive constant c, independent of Z. On the other hand, also using this simple 
scaling relation one sees that 

K(^) = dZ^ 3 , 

for some positive constant d (independent of Z). Thus, again we observe the same dependence 
in the number of particles as in the previous example, i.e., 



TF 



z 1 / 3 ' 

Thus, for large values of Z, the second and third terms in fljjD are negligible. 



8 
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